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Precision measurement of 



the 




oscillation frequency with the decay 



A key ingredient to searches for physics beyond the Standard Model in B® mixing 
phenomena is the measurement of the B®-B® oscillation frequency, which is equiv- 
alent to the mass difference Am s of the B® mass eigenstates. Using the world's 
largest B® meson sample accumulated in a dataset, corresponding to an integrated 
luminosity of 1.0 fb _1 , collected by the LHCb experiment at the CERN LHC in 2011, 
a measurement of Am s is presented. A total of about 34,000 B® — > D~ir + signal 
decays are reconstructed, with an average decay time resolution of 44 fs. The oscilla- 
tion frequency is measured to be Am s = 17.768 ± 0.023 (stat) ± 0.006 (syst) ps _1 , 
which is the most precise measurement to date. 
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1 Introduction 



The Standard Model (SM) of particle physics, despite its great success in describing 
experimental data, is considered an effective theory only valid at low energies, below the 
TeV scale. At higher energies, new physics phenomena are predicted to emerge. For 
analyses looking for physics beyond the SM (BSM) there are two conceptually different 
approaches: direct and indirect searches. Direct searches are performed at the highest 
available energies and aim at producing and detecting new heavy particles. Indirect 
searches focus on precision measurements of quantum-loop induced processes. Accurate 
theoretical predictions are available for the heavy quark sector in the SM. It is therefore 
an excellent place to search for new phenomena [lj[2], since any deviation from these 
predictions can be attributed to contributions from BSM. 

In the SM, transitions between quark families (flavours) are possible via the charged 
current weak interaction. Flavour changing neutral currents (FCNC) are forbidden at 
lowest order, but are allowed in higher order processes. Since new particles can contribute 
to these loop diagrams, such processes are highly sensitive to contributions from BSM. An 
example FCNC transition is neutral meson mixing, where neutral mesons can transform 
into their antiparticles. Particle-antiparticle oscillations have been observed in the K°-K° 
system (3), the B°-B°_Q system, the B°-B° s system |H|[6) and the D°-D° system (7 



fO 



The frequency of B®-B® oscillations is the highest. On average, a B® meson changes its 
flavour nine times between production and decay. This poses a challenge to the detector 
for the measurement of the decay time. Another key ingredient of this measurement is the 
determination of the flavour of the B° s meson at production, which relies heavily on good 
particle identification and the separation of tracks from the primary interaction point. 

The observed particle and antiparticle states B® and B° s are linear combinations of 
the mass eigenstates B H and B L with masses m H and and decay widths T H and 



Tl, respectively 11 . The B® oscillation frequency is equivalent to the mass difference 
Am s = me — ^l- The parameter Am s is an essential ingredient for all studies of time- 
dependent matter-antimatter asymmetries involving B° s mesons, such as the B° s mixing 
phase cj) s in the decay B® — > J/ip<f) [l2j. It was first observed by the CDF experiment j(3 j. The 
LHCb experiment published a measurement of this frequency using a dataset, corresponding 



to an integrated luminosity of 37 pb -1 , taken in 2010 13 . This analysis complements the 
previous result and is obtained in a similar way, using a data sample, corresponding to an 
integrated luminosity of 1.0 fb~ , collected by LHCb in 2011. 



2 The LHCb experiment 

The LHCb experiment is designed for precision measurements in the beauty and charm 
hadron systems. At a center-of-mass energy of \fs = 7 TeV, about 3 ■ 10 11 bb pairs were 
produced in 2011. The LHCb detector [14] is a single-arm forward spectrometer covering 
the pseudorapidity range from two to five. The excellent decay time resolution necessary to 
resolve the fast B° s -B° s oscillation is provided by a silicon-strip vertex detector surrounding 
the pp interaction region. At nominal position the sensitive region of the vertex detector 
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is only 8 mm away from the beam. Impact parameter (IP) resolution of 20 urn for tracks 
with high transverse momentum (pp) is achieved. 

Charged particle momenta are measured with the LHCb tracking system consisting of 
the aforementioned vertex dector, a large-area silicon-strip detector located upstream of 
a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes placed downstream. The combined tracking system has 
momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c. 

Since this analysis is performed with decays involving only hadrons in the final state, 
excellent particle identification is crucial to suppress background. Charged hadrons 



are identified using two ring-imaging Cherenkov detectors 15 . Photon, electron, and 
hadron candidates are identified by a calorimeter system consisting of scintillating-pad and 
preshower detectors, an electromagnetic calorimeter, and a hadronic calorimeter. Muons 
are identified by a system composed of alternating layers of iron and multiwire proportional 
chambers. 

The first stage of the trigger 16 is implemented in hardware, based on information from 
the calorimeter and muon systems, and selects events that contain candidates with large 
transverse energy and transverse momentum. This is followed by a software stage which 
applies a full event reconstruction. The software trigger used in this analysis requires a 
two-, three- or four-track secondary vertex with a significant displacement from the primary 
interaction, a large sum of pr of the tracks, and at least one track with pp > 1.7 GeV/c. In 
addition an IP x 2 with respect to the primary interaction greater than 16 and a track fit 
X 2 per degree of freedom < 2 is required. The IP x 2 is defined as the difference between 
the x 2 °f the primary vertex reconstructed with and without the considered track. A 
multivariate algorithm is used for the identification of the secondary vertices. 

For the simulation, pp collisions are generated using Pythia 6.4 [17] with a specific 
LHCb configuration [18] . Decays of hadronic particles are described by EvtGen fl9] , 



in which final state radiation is generated using Photos 20 . The interaction of the 



generated particles with the detector and its response are implemented using the Geant4 



toolkit 21 , as described in Ref. 



22 



3 Signal selection and analysis strategy 

The analysis uses B® candidates reconstructed in the flavour-specific decay moda^ 
B° s ^ DJ7T+ in five D~ decay modes, namely D~ ->• <p(K + K-)n- , D~ — > K*°(K + n')K~ , 
D~ — > K + K~tt~ nonresonant, D~ — > K~7t + 7t~, and D~ — > 7t~7t + tt~. To avoid double 
counting, events that contain a candidate passing the selection criteria of one mode, are 
not considered for the subsequent modes, using the order listed above. All reconstructed 
decays are flavour-specific final states, thus the flavour of the B® candidate at the time 
of its decay is given by the charges of the final state particles. A combination of tagging 
algorithms is used to identify the B® flavour at production. The algorithms provide for 

1 Unless explicitly stated, inclusion of charge-conjugated modes is implied. 
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each candidate a tagging decision as well as an estimate of the probability that this decision 
is wrong (mistag probability) . These algorithms have been optimized using large event 



samples of flavour-specific decays 23 , 24 



To be able to study the effect of selection criteria that influence the decay time 
spectrum, we restrict the analysis to those events in which the signal candidate passed 
the requirements of the software trigger algorithm used in this analysis. Specific features, 
such as the masses of the intermediate and K*° resonances or the Dalitz structure of the 
D~ — > 7t~7t + tt~ decay mode, are exploited for the five decay modes. The most powerful 
quantity to separate signal from background common to all decay modes is the output 



of a boosted decision tree (BDT) |25 . The BDT exploits the long B° s lifetime by using 



as input the IP x °f f ne daughter tracks, the angle of the reconstructed B s momentum 
relative to the line between the reconstructed primary vertex, and the B° s vertex and the 
radial flight distance in the transverse plane of both the B® and the D~ meson. Additional 
requirements are applied on the sum of the pt of the B° s candidate's decay products as well 
as on particle identification variables, and on track and vertex quality. The reconstructed 



D s mass is required to be consistent with the known value 26 . After this selection, a 
total of about 47,800 candidates remain in the B® — > D~n + invariant mass window of 
5.32 - 5.98 GeV/c 2 . 

An unbinned likelihood method is employed to simultaneously fit the B® invariant mass 
and decay time distributions of the five decay modes. The probability density functions 
(PDFs) for signal and background in each of the five modes can be written as 



V = V m (m) V t (t, q\a t , rj) V at (a t ) V v {n), (1) 

where m is the reconstructed invariant mass of the B° s candidate, t is its reconstructed 
decay time and a% is an event-by-event estimate of the decay time resolution. The tagging 
decision q can be if no tag is found, —1 for events with different flavour at production and 
decay (mixed) or +1 for events with the same flavour at production and decay (unmixed). 
The predicted event-by-event mistag probability rj can take values between and 0.5. 
The functions V m and Vt describe the invariant mass and the decay time probability 
distributions, respectively. Vt is a conditional probability depending on a t and r\. The 
functions V at and V v are required to ensure the proper relative normalization of Vt for 
signal and background [27]. The functions V at and V v are determined from data, using 
the measured distribution in the upper B® invariant mass sideband for the background 
PDF and the sideband subtracted distribution in the invariant mass signal region for the 
signal PDF. 

This measurement has been performed "blinded" , meaning that during the analysis 
process the fitted value of Am s was shifted by an unknown value, which was removed 
after the analysis procedure had been finalized. 
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Figure 1: Invariant mass distributions for D~ir + candidates with the D~ meson de- 

caying as a) D~ -> 4>(K + K~)tt~ , b) Dj ^ K*°(K + ir~)K~ , c) D~ -> K + K~tt~ nonresonant, 
d) D J— > K~7r + ir~, and e) D~ — > 7r~7r + 7r~. The fits and the various background components 
are described in the text. Misidentified backgrounds refer to background from B° and A9 decays 
with one misidentified daughter particle. 



4 Invariant mass description 

The invariant mass of each candidate is determined in a vertex fit constraining the 



s 

D~ invariant mass to its known value 26 . The invariant mass spectra for the five decay 
modes after all the selection criteria are applied are shown in Fig. [TJ The fit to the five 
distributions takes into account contributions from signal, combinatorial background and 
fe-hadron decay backgrounds. The signal components are described by the sum of two 
Crystal Ball (CB) functions [28], which are constrained to have the same peak parameter. 
The parameters of the CB function describing the tails are fixed to values obtained 
from simulation, whereas the mean and the two widths are allowed to vary. These are 
constrained to be the same for all five decay modes. It has been checked on data that the 
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Table 1: Number of candidates and B® signal fractions in the mass range 5.32 - 5.98 GeV/c 2 . 



Decay mode (D s vr + ) candidates f B o^ D - n + fB°^-D T K ± 



D~ ^ (f)(K + K-)tt- 


14691 


0.834 ± 0.008 


D~ — >■ K*°(K + tt~)K~ 


10866 


0.857 ± 0.009 


D~ — > K + K~ir~ nonresonant 


11262 


0.595 ± 0.009 


D~ -> K~n + n~ 


4288 


0.437 ± 0.014 


D~ — > ir~ir + iT~ 


6674 


0.599 ± 0.008 0.019 ±0.010 



Total 47781 0.714 ± 0.004 0.019 ±0.010 



mass resolution is compatible among all modes. 

The 6-hadron decay background includes B° and A® decays with one misidentified 
daughter particle. Their mass shapes are derived from simulated samples. The yields for 
the different 6-hadron decay backgrounds are allowed to vary individually for each of the 
five decay modes. Another component originates from — > DJK^ 1 decays, in which the 
kaon is misidentified as a pion. This contribution is treated as signal in the decay time 
analysis. 

The requirement that the invariant mass be larger than 5.32 GeV/c 2 rejects background 
candidates from B° s decays with additional particles in the decay not reconstructed, such 
as -B°— >■ D*~tc + (D*~ — > D~7i° or D~ 7). The fitted number of signal candidates does not 
change with respect to a fit in a larger mass window. The high mass sideband region 
5.55 - 5.98 GeV/c 2 provides a sample of mainly combinatorial background candidates. 
The mass distribution is described by an exponential function, whose parameters are 
allowed to vary individually for the five decay modes. By including this region in the fit, 
we are able to determine the decay time distribution as well as the tagging behaviour of 
the combinatorial background. 

The number of used candidates, along with the signal fractions extracted from the 
two dimensional fit in mass and decay time, are reported in Table [TJ One complication 
arises from the fact that the shape of the invariant mass distribution of the B° s Df^ 
events is very similar to that of the B° background. Therefore the fraction of B® — > DJK^ 1 
candidates has been determined in a fit to the D~ — > 7t~7t + tt~ mode only, in which no B° 
background is present. Subsequently this value is used for all other modes. 



5 Decay time description 

The decay time of a particle is measured as 

Lm , , 

t = — , 2 
V 

where L is the distance between the production vertex and the decay vertex of the particle, 
m its reconstructed invariant mass, and p its reconstructed momentum. We use the decay 
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time calculated without the D~ mass constraint to avoid a systematic dependence of the 
Bg decay time on the reconstructed invariant mass. The theoretical distribution of the 
decay time, t, ignoring the oscillation and any detector resolution, is 

V t oc T s e- Tat cosh f -^t j 6(t), (3) 



where F s is the B® decay width and Ar s the decay width difference between the light 
and heavy mass eigenstatej^] The value for Ar s is fixed to the latest value measured 



by LHCb [12] Ar s = 0.106 ± 0.011 ± 0.007 ps" 1 . It is varied within its uncertainties to 
assess the systematic effect on the measurement of Am s . The Heaviside step function 8(t) 
restricts the PDF to positive decay times. 

To account for detector resolution effects, the decay time PDF is convolved with a 
Gaussian distribution. The width a t is taken from an event-by-event estimate returned 
by the fitting algorithm that reconstructs the B° s decay vertex. Due to tracking detector 
resolution effects at needs to be calibrated. A data-driven method, combining prompt D~ 
mesons from the primary interaction with random ir + mesons, forms fake B® candidates. 
The decay time distribution of these candidates, each divided by its event-by-event at, is 
fitted with a Gaussian function. The width provides a scale factor S at = 1.37, by which 
each a t is multiplied, such that it represents the correct resolution. By inspecting different 
regions of phase space of the fake B° s candidates, the uncertainty range on this number is 
found to be 1.25 < S at < 1.45. The variation is taken into account as part of the Am s 
systematic studies. The resulting average decay time resolution is S at x (a t ) = 44 fs. 

Some of the selection criteria influence the shape of the decay time distribution, e.g. the 
requirement of a large IP for B° s daughter tracks. Thus a decay time acceptance function 
St(t) has to be taken into account. Its parametrization is determined from simulated data 
and the parameter describing its shape is allowed to vary in the fit to the data, while T s is 
fixed to the nominal value [26|. Taking into account resolution and decay time acceptance, 



the PDF given in Eq.([3]) is modified to 



Vt(t\cr t ) oc 



r s e" rsi cosh i^t ) 6(t) 



G(t;0,S at a t ) 8 t (t), (4) 



V 2 

with G(t; 0, S at a t ) being the resolution function determined by the method mentioned 
above. The decay time PDFs for the B° and A® backgrounds are identical to the signal 
PDF, except for Ar being zero, and T s being replaced by their respective decay widths (26]. 
The shape of the decay time distribution of the combinatorial background is determined 
with high mass sideband data. It is parametrized by the sum of two exponential functions 
multiplied by a second order polynomial distribution. The exponential and polynomial 
parameters are allowed to vary in the fit and are constrained to be the same for the five 
decay modes. 



2 AT S and Am s are measured in units with h— 1 throughout this paper. 
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6 Flavour tagging 



To determine the flavour of the B® meson at production, both opposite-side (OST) and 
same-side (SST) tagging algorithms are used. The OST exploits the fact that b quarks at 
the LHC are predominantly produced in quark-ant iquark pairs. By partially reconstructing 
the second b hadron in the event, conclusions on the flavour at production of the signal B® 
candidate can be drawn. The OST have been optimized on large samples of B + — > J/tpK + , 
B -> n + D*~X, and B° ->■ D-tt + decays (23) . 

The SST takes advantage of the fact that the net strangeness of the pp collision is 
zero. Therefore, the s quark needed for the hadronization of the B® meson must have been 
produced in association with an s quark, which in about 50% of the cases hadronizes to 
form a charged kaon. By identifying this kaon, the flavour at production of the signal B° s 
candidate is determined. The optimization of the SST was performed on a data sample of 
B®—> D~7r + decays, which has a large overlap with the sample used in this analysis 24 



However, since the oscillation frequency is not correlated with the parameters describing 
tagging performance, this does not bias the Am s measurement. 

The decisions given by both tagging algorithms have a probability u to be incorrect. 
Each tagging algorithm provides an estimate for the mistag probability rj which is the 
output of a neural network combining various event properties. The true mistag probability 



u> can be parametrized as a linear function of the estimate rj 23 , 24 



u)=p +pi x (77- (77)), (5) 

with (77) being the mean of the distribution of 77. This parametrization is chosen to 
minimize the correlations between p Q and p\. The calibration is performed separately for 
the OST and SST. 

The sets of calibration parameters (po>Pi)ost and (po>.Pi)sst are allowed to vary in the 
fit. The figure of merit of these tagging algorithms is called the effective tagging efficiency 
£ e g. It gives the factor by which the statistical power of the sample is reduced due to 
imperfect tagging decisions. In this analysis, e e g is found to be (2.6 ± 0.4)% for the OST 
and (1.2 ± 0.3)% for the SST. Uncertainties are statistical only. 



7 Measurement of Am. 



Adding the information of the flavour tagging algorithms, the decay time PDF for tagged 
signal candidates is modified to 



V t {t\a t ) oc \Y s e 



-T s t 



cosh 



t) + q [1 - 2w(t7ost, Vsst)] cos(Am s t) 



9(t) 



G(t,S at cr t ) S t (t) e, 



(6) 



where e gives the fraction of candidates with a tagging decision. Signal candidates without 
a tagging decision are still described by Eq.Q multiplied by an additional factor (1 — e) 
to ensure the relative normalization. 
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Figure 2: Decay time distribution for the sum of the five decay modes for candidates tagged as 
mixed (different flavour at decay and production; red, continuous line) or unmixed (same flavour 
at decay and production; blue, dotted line). The data and the fit projections are plotted in a 
signal window around the reconstructed B® mass of 5.32 - 5.55 GeV/c 2 . 



The information provided by the opposite-side and same-side taggers for the signal is 
combined to a single tagging decision q and a single mistag probability w(r]osT, Vsst) using 
their respective calibration parameters Po OS t/sst an< ^ Piost/sst- The individual background 
components show different tagging characteristics for candidates tagged by the OST or 
SST. The b hadron backgrounds show the same opposite-side tagging behaviour (q and 
u) as the signal, while the combinatorial background shows random tagging behaviour. 
For same-side tagged events, we assume random tagging behaviour for all background 
components. We introduce tagging asymmetry parameters to allow for different numbers 
of candidates being tagged as mixed or unmixed, and other parameters to describe the 
tagging efficiencies for these backgrounds. As expected, the fitted values of these asymmetry 
parameters are consistent with zero within uncertainties. 

All tagging parameters, as well as the value for Am s , are constrained to be the same 
for the five decay modes. The result is Am s = 17.768 ± 0.023 ps^ 1 (statistical uncertainty 
only). The likelihood profile was examined and found to have a Gaussian shape up to 
nine standard deviations. The decay time distributions for candidates tagged as mixed 
or unmixed are shown in Fig. [2j together with the decay time projections of the PDF 
distributions resulting from the fit. 
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Table 2: Systematic uncertainties on the Am s measurement. The total systematic uncertainty 
is calculated as the quadratic sum of the individual contributions. 

Source Uncertainty [ps -1 ] 



z-scale 0.004 

Momentum scale 0.004 

Decay time bias 0.001 

Total systematic uncertainty 0.006 



8 Systematic uncertainties 



With respect to the first measurement of Am s at LHCb 13 all sources of systematic 
uncertainties have been reevaluated. 

The dominant source is related to the knowledge of the absolute value of the decay time. 
This has two main contributions. First, the imperfect knowledge of the longitudinal (z) 
scale of the detector contributes to the systematic uncertainty. It is obtained by comparing 
the track-based alignment and survey data and evaluating the track distribution in the 
vertex detector. This results in 0.02% uncertainty on the decay time scale and thus an 
absolute uncertainty of ±0.004 ps -1 on Am s . 

The second contribution to the uncertainty of the decay time scale comes from the 
knowledge of the overall momentum scale. This has been evaluated by an independent 
study using mass measurements of well-known resonances. Deviations from the reference 



values 26 are measured to be within 0.15%. However, since both the measured invariant 
mass and momentum enter the calculation of the decay time, this effect cancels to some 
extent. The resulting systematic on the decay time scale is evaluated from simulation to 
be 0.02%. This again translates to an absolute uncertainty of ±0.004 ps -1 on Am s . 

The next largest systematic uncertainty is due to a possible bias of the measured decay 
time given by the track reconstruction and the selection procedure. This is estimated 
from simulated data to be less than about 0.2 fs, and results in ±0.001 ps _1 systematic 
uncertainty on Am s . 

Various other sources contributing to the systematic uncertainty have been studied 
such as the decay time acceptance, decay time resolution, variations of the value of Ar s , 
different signal models for the invariant mass and the decay time resolution, variations 
of the signal fraction and the fraction of B® — > DJK^ 1 candidates. They are all found 
to be negligible. The sources of systematic uncertainty on the measurement of Am s are 
summarized in Table [2] 

9 Conclusion 

A measurement of the B° s -B° s oscillation frequency Am s is performed using B° s — > Djtt + 
decays in five different D~ decay channels. Using a data sample corresponding to an 
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integrated luminosity of 1.0 fb 1 collected by LHCb in 2011, the oscillation frequency is 
found to be 



Am s = 17.768 ± 0.023 (stat) ± 0.006 (syst) ps~\ 
in good agreement with the first result reported by the LHCb experiment 



current world average, 17.69 ± 0.08 ps 26 . This is the most precise measurement of 



13 and the 



Am s to date, and will be a crucial ingredient in future searches for BSM physics in B® 
oscillations. 
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